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We present the results of Monte Carlo simulations of two different 
10-state Potts glasses with random nearest neighbor interactions on 
a simple cubic lattice. In the first model the interactions come from 
a ± J distribution and in the second model from a Gaussian one, and 
in both cases the first two moments of the distribution are chosen 
to be equal to Jo = —1 and A J = 1. At low temperatures the spin 
autocorrelation function for the ±J model relaxes in several steps 
whereas the one for the Gaussian model shows only one. In both 
systems the relaxation time increases like an Arrhenius law. Unlike 
the infinite range model, there are only very weak finite size effects 
and there is no evidence that a dynamical or a static transition exists 
at a finite temperature. 

PACS numbers: 64.70.Pf, 75.10.Nr, 75.50.Lk 

In recent years generalized spin-glass-type models such as the p-spin model with 
p > 3 or the Potts glass with p > 4, have found a large attention |2|, [|, ^ 
||, [7|, H ||, |10|, [11], [T|, [13| as prototype models for the structural glass transition 



[15], [IB, [T7|, |T^, [19]. In the case of infinite range interactions, i.e. mean 
field, these models can be solved exactly and it has been shown that they have a 
dynamical as well as a static transition at a temperature To and To, respectively 
0, |, |, |, |, 1, 0, i i H [J2 H [U. At T D the relaxation times diverge but 



no singularity of any kind occurs in the static properties, whereas at To a nonzero 
static glass order parameter appears discontinuously. Close to To the time and 
temperature dependence of the spin autocorrelation function is described by the 
same type of mode coupling equations || 0] that have been proposed by the 
idealized version of mode coupling theory (MCT) |15|, [16j for the structural glass 



transition which suggests a fundamental connection between these rather abstract 
spin models and real structural glasses. 



Now it is well known that for real glasses the divergence of the relaxation times 
predicted by MCT is rounded off since thermally activated processes, which are not 
taken into account by this version of the theory, become important |HJ. This 
is in contrast to the mean field case because there these processes are completely 
suppressed since the barriers in the (free) energy landscape become infinitely high 
if T < Tjj. To what extend the static transition that exists in mean field can be 
seen also in the real system, is a problem whose answer is still controversial 17, 



m, mi 



In the present paper, we investigate whether the transitions present in the case 
of a mean field 10-state Potts glass are also found if the interactions are short 
ranged. The mean field model has recently been studied in great detail by means 
of Monte Carlo methods and the results are in good qualitative agreement with 
the main predictions of the one-step replica symmetry breaking theory fl], 
although very strong finite size effects occur, the nature of which are not fully 



understood [PL2| . 



The Hamiltonian of the short range model is 

w = -5>«(p<W,-i)- (i) 

The Potts spins <7j on the lattice sites i of the simple cubic lattice take p discrete 
values, Oi G {1,2, .. . ,p}, and the index j runs over the six nearest neighbors of 
site i. The exchange constants Jjj are taken either from a bimodal distribution 

P(J l3 ) = xSiJij - J) + (1 - x)5(J ij + J) (2) 

or a Gaussian distribution 

P(J ii ) = [VMA J)]- 1 exp{-(J ii - J ) 2 /2(A J) 2 }. (3) 

In both cases the first two moments are chosen to be Jo = [Jij]av = — 1 and 
(A J) 2 = [Jfj]av — [Jij]tv = +1- (F° r P = 10, a sufficiently negative J is necessary to 
avoid ferromagnetic order and we have indeed found that neither the magnetization 
nor the magnetic susceptibility show any sign of ferromagnetic ordering ||li|).) F° r 
the distribution given by Eq. (fj) this choice means J = V2 and x = (l-l/V2)/2 w 
0.146. We carry out Monte Carlo runs with the standard heat-bath algorithm fl2"0| 
making up to 10 8 Monte Carlo steps (MCS) per spin (for equilibration as well 
as production), with lattices of linear dimensions L = 6, L = 10 and L = 16, 
and using periodic boundary conditions. The average over the quenched disorder 



[• • • ] av is realized by averaging over 100 independent realizations of the system for 
L = 6 and 10 and over 50 realizations for L = 16. For T = 1.3 and 1.4 and L = 10 
only 10 realizations were used. In the following we will set the Boltzmann constant 
&b = 1 and measure temperature in units of AJ. 

In Fig. m we show the temperature dependence of the energy per spin e(T) and 
of the specific heat c(T), for the model given by Eq. (0). Both quantities seem to 
be essentially independent of system size, in stark contrast to the results for the 
mean field case ||12|| . Furthermore we see that for T < 2 the energy is basically 
constant and we have found that at low T the temperature dependence of the 
specific heat is of the form (A/T) 2 exp(— A/T), with A « 14, i.e. a dependence 
expected for a two-level system with asymmetry A. Thus we conclude that this 
data does not show any evidence that in the temperature range investigated a 
static transition occurs. 

A similar conclusion is reached from the T— dependence of the spin glass sus- 
ceptibility xsg — N[{q 2 )] av / (j) — 1). Here the symmetrized order parameter q is 
defined as (11], [12] 



\ 



1 p_1 

-J2(^ U ) 2 with q^^L-zJ^iS^nS^Y, n,v = l,2,... ,p-l 

p — 1 z — 4 z — ' 

(4) 
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and the spins Si a refer to the simplex representation of the Potts spin <7j in replica 
a [^] . (Note that in order to define an order parameter it is necessary to consider 
two replicas of the system, i.e. two systems with the same realization of the 
disorder.) Figure |]a shows the T— dependence of xsg an d we see that this quantity 
is almost independent of L. Furthermore we recognize that xsg remains finite 
even as T — > 0. Figure 0b shows the (scaled) first moment of q and we conclude 
that it decreases like 1/a/Z^, i.e. shows a trivial size dependence for all T. Also 
an analysis of fourth order cumulants (not shown here) confirms the conclusion 



that no static transition occurs |]T3|. Similar results are found for the model with 
Gaussian distribution, Eq. @. 

In order to see whether this model has a dynamic transition at a finite temper- 
ature we consider the time autocorrelation function of the Potts spins: 

C(t) = [L 3 (p-l)]- 1 Y,[(S i (t')-S l (t' + t))} av . (5) 

i=l 

The time dependence of C(t) is shown in Fig. |||a for various temperatures. We see 
that at intermediate temperatures C{t) has a plateau at around 0.6 with a width 



that increases quickly with decreasing T. Such a behavior is very reminiscent of 
the time and T dependence of glass forming systems close to the MCT temperature 
Td B 111 0, [13|, 14, 15, IT]. However, as we will show below, in this case the 



reason for the existence of a plateau is very different. Interestingly enough C(t) 
shows at low temperatures a second plateau, and also its length increases rapidly 
with decreasing temperature. Such a multi-step relaxation has so far been seen 
only for few glass forming systems |E| and is a rather unusual behavior. Below we 
will come back to this feature and discuss its origin in more detail. That the time 
dependence of C(t) is basically independent of the system size is demonstrated in 
the inset of Fig. |3]a. Also this behavior differs strongly from the one found for the 
mean field model ]12| . 



In order to study the slowing down of the dynamics of the system we can 
define relaxation times r« via C(ri) = 0.4, Cfo) = 0.08 and C{t^) = 0.05. These 
definitions characterize the relaxation times for the different processes seen in 
Fig. |3|a. The temperature dependence of Tj is shown in Fig. ||b. We see that this 
dependence can be described very well by an Arrhenius law with an activation 
energy around 28.2 and 14.6 (straight lines). Such a T— dependence is not in 
agreement with the expectation from MCT which predicts around a power-law 
dependence. The observation that the activation energies are close to p\J\ ~ 14.1 
and 2p\J\ ~ 28.3 suggests instead that the relaxation of the spins is given by the 
breaking of one and two bonds, respectively. This interpretation is also supported 
by the time dependence of the autocorrelation function of individual spins ||13|| , 
since these functions typically fall in three classes: Those that are relaxing fast 
(spins that have only negative bonds), those that relax on intermediate times (one 
bond needs to be broken) and those that relax slowly (two bonds that have to be 
broken). Using these arguments and the concentration of ferromagnetic bonds, x, 
one can estimate also the height of the plateaus which is predicted to be 0.61 



and 0.13, which is in very good agreement with the height that can be estimated 
from Fig. |3]a to be 0.59 and 0.12, respectively. 

We emphasize that the occurrence of several plateaus in the spin autocorre- 
lation function (Fig. |3]a) is an immediate consequence of the bond distribution, 
Eq. (H), and hence no such plateaus are expected to occur for the Gaussian bond 
distribution, Eq. (|3]). This expectation is indeed born out by the numerical data 
shown in Fig. |]a since C(t) does not show any sign for a plateau for any of the 
temperatures investigated. If one uses again various definitions of relaxation times 



Tj(T), this time the values 0.1, 0.05 and 0.02, their T— dependence is found again 
to be Arrhenius like, see Fig. |]b. However, in contrast to the results for the ±J 
distribution the activation energies depend continuously on the definition of r^. 
Thus one finds that the dynamical transition that the same model exhibits for the 
infinite range of interactions is completely wiped out for the nearest neighbor case 
and we conclude that neither a remnant of the static glass transition nor of the 
dynamic transition exists. 

Obviously the dynamical behavior of the present model is not similar to the one 
of a supercooled liquid close to its glass transition. But such a similarity can be 
expected to occur if we consider a variant of the model that interpolates between 
the short range version of the Potts glass and its infinite range version. E.g. we 
could choose interactions that have a finite range but that extend further than the 
nearest neighbors. For such a model we expect that a "rounded" version of the 
dynamic transition will occur, i.e., the autocorrelation function develops a long- 
lived plateau and the relaxation time exhibits the onset of a power law singularity, 
before it crosses over to a simple Arrhenius behavior. This is the behavior found 
for large but finite mean field systems @, 0, HI El- Much less 

can be expected 



to be seen in the static properties of that model: For any finite range of the 
interaction, no nonzero order parameter can appear H, and a jump singularity 



in xsg cannot occur either. All what remains is a rounded kink in the entropy vs. 
temperature curve at the rounded static transition, thus avoiding the Kauzmann 
[24|j catastrophe of a vanishing configurational entropy. In fact, this description is 
nicely consistent with all known facts about the structural glass transition. Thus, 
if the analogy between the latter and the behavior of medium range Potts glasses 
goes through, the search for a static glass transition will remain elusive! This is 
a somewhat surprising conclusion, since for p = 2 the model reduces to the Ising 
spin glass, where one knows that a (second-order) glass transition temperature 



occurs at T c > [25|, and also for the p = 3 Potts glass there seems to be at least 



a divergent xsg as T —>■ |26|, p7| . However, for large values of p the static and 
dynamic properties seem to be very different. 
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FIG. 1: Temperature dependence of the energy per spin e(T), (a), and specific heat per 
spin, (b). The different curves correspond to different lattice sizes L. 
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FIG. 2: Temperature dependence of the spin glass susceptibility XSG-, (a)> an d scaled 
first moment of the order parameter [(q)] a v, (b). The different curves correspond to 
different lattice sizes L. 




FIG. 3: (a) Time dependence of the spin autocorrelation function C(t) for temperatures 
(left to right) T = 2.4, 2.2, 2.0, 1.8, 1.7, 1.6, 1.5, 1.4, and 1.3 (L = 10). The horizontal 
dashed lines are used to define the relaxation times Tj (see text for details). In the inset 
we compare C(t) for L = 10 and L = 16 at T = 1.6. (b) Symbols: Arrhenius plot of 
the relaxation time Tj (see text for definition). Straight lines: Fits to the data with an 
Arrhenius law. The numbers are the activation energies. 
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FIG. 4: (a) Time dependence of the spin autocorrelation function C(t) for temperatures 
T = 2.8, 2.6, 2.4, 2.2, and 2.0 (left to right) for the Gaussian model. The horizontal 
dashed lines are used to define the relaxation times Tj (see text for details), (b) Symbols: 
Arrhenius plot of the relaxation time Ti (see text for definition). Straight lines: Fits to 
the data with an Arrhenius law. The numbers are the activation energies. 



